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Although many tables of integrals occurring in quantum chemistry work have 
been published [1, 2, 3], their application is, as a rule, limited to calculations 
involving Slater functions, so that the use of Hartree-Fock self-consistent field 
(SCF) functions in analytical form necessitates recourse to other means of com- 
putation. The early Slater method [4] of transforming numerical SCF functions 
into analytical form has been greatly improved by Léwdin & Appel [5], who 
introduced a method of successive approximations for the expansion of a given 
SCF function f(n,1)=r-R(n,l) in exponentials 


f(n, Dax I (r;—1) 2 An exp (ant); [ro )- dra Fis 
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where R(n,1)=the radial part of the complete wave function; 
n=the principal quantum number; 
1=the azimutal quantum number; 
;=index taking the integral values satisfying the inequality 
0<i<n-—I—-1; 

7,=the ith zero point of f(n, 1); 

Am, 4m, = constants giving the best fit; 

C (n, l) =normalization constant. 


Making use of the Aitken iteration scheme [6] for a simple and accurate 
determination of the zero points 7; of the numerical SCF functions, the tabu- 
lated values, after division by the polynomial r'*1 TT], (4, —1), represent a node- 
less and monotonous function, expressible as a sum of exponentials. The method, 
correctly used, leads to an almost perfect fit (fanai — ftavw = +£0.001-002). It can 
be shown that these analytical wave functions give rise to the same general 
form of integrals as those involved in calculations with Slater functions, the 
contributory terms, however, differing in magnitude in the two cases. 

In this note an electronic computer (BESK) program is described, which per- 
mits the computation of the following one-electron two centre integrals 
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Overlap integral: [vs wedt Natnp sD 


Moment integral: | Wazppat 
; 1 1 
Resonance integral: | vs cyedes | vacvade 
TA 'p 


Kinetic energy integral: | wa(—4A)yppdt 


1 1 
Nuclear attraction integral: | Yaa dt; | pe YB dt, 
B A 


where the indices refer to atoms A and B respectively, according to Fig. 


1. 


The rather modest restriction imposed on the principal quantum numbers n, 


and ng can, if needed, be removed by simple changes of the program. 
The calculations are as usual made by means of spheroidal coordinates 
atomic units having the following properties 


1<&<oco o=internuclear distance 


¢=the azimutal angle 


fata 40 (en) 

r4 cos 04=4$0 (1+ én) 

rp COS p= to (1—En) 

ra sin 04=1g sin O2=40V(E—1) (77-1) 
dt=to° (&—n*)dédndd. 


in 


The various molecular integrals considered in this work can be given the com- 


mon form 
Jy=Ky ky, ye Ch [Ay (Vx) wee (6,)]a, 


where N= a number labelling the integral 


Ky, ky, C2, Yu, 6x = constants 


Ax(yx)= | "exp (—yx€) dé, 


1 
B, (6x) = eaters (—6,) dn 
1 
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Fig. 1. 


and [A,,(y,):B,(d,)] represents a sum of products A, B,, the rule of summation 
being determined by the complete expression for the particular integral under 
consideration (cf. example below). 

The well-known molecular functions A, (y) and B,(6) are given in tables of 
molecular integrals and may be calculated in different ways (see e.g. Kotani 
et al. [1], Corbatdé [7], Flodmark [8], Miller et al. [3]). In calculating the B, (6) 
function a difficulty arises due to the slow convergence of this function for low 
values of |6/y|. This has been avoided in the present work by the use of a 
recursion formula for high values of |5/y| and a series expansion for low ones. 
Details are given in Appendix I. 

In order to exemplify an explicit form of the common integral expression we 
shall consider one of the integrals, which appear in calculating the overlap 
between hybrid orbitals of hydrogen and fluorine. Thus consideration of the 
1s-orbital of hydrogen and the 2s-orbital of fluorine gives the following integral 


Sy= | vats) ye@s)ar. 


The complete wave functions are expressed by 


= 


va(ls)= =e 


{-1 
22 pn UT | fr (28); 


fx (2s) being a normalized SCF function given analytical form according to 
Léwdin & Appel [5] from numerical data given by Brown [9]. 
fe (28) = Cp (28) +7 (0.25358 — 1) {3.6103 exp (— 1.9323 r) + 22.566 - exp (— 3.4577 r) 
+ 18.600-exp (—6.43367)} 


3 
= Op (28)-r (7-7) X Am exp (—a,°7r). 
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Arbitrarily putting the hydrogen atom as the A-atom we get 


1 3 
Sy = Cz (28) ss | ex alee (rin — 1p) 2, An exp (—Qm°%z) OT; 


which is transformed to spheroidal coordinates with the aid of the relations 
r= 2 (E+) te =S(E-n) dr =F (Pn?) dE dn dg, 


where 9 is the internuclear distance of the hydrogen-fluorine bond under con- 
sideration. 
Foetus) 


m=1 


3 
2 A {rip (A, (Ym) r Bo (Om) * Ay (Ym) . B, (Om)) 
— 40 ((Ag (Ym) * Bo (Om) + Ao (Ym) *-Bs (Sm) — Ay (Ym) * Bs (8m) — Ag (7m) Bs (Om))} 


3 2 
= Ky 33 An PJ a [Ay (1m) - B, (Om) Ja 


m=1 


which is of the form given above. The parameter tape will contain the values 
of N, Ky, k,,c, uw and y, y, and 6, together with some parameters giving the 
sign of the terms in the integral. 

Using the expression for the overlap, it is easy to write down the corre- 
sponding moment integral. As shown in appendix II the moment integral My 
is related to the overlap integral Sy(A,(y)-B,(6)) in the following way 


My = $ @ Sy(Aps1 (y) + By +1 (0). 


Because of this very simple relationship the program has been extended auto- 
matically to perform the computation of the moment integral from the para- 
meter values given for the overlap integral. In a program used for the calcula- 
tion of molecular dipole moments (to be published), the procedure saved half 
of the stamping work. The computations of overlap and moment integrals in- 
volving s and p functions according to this program have been carried out on 
the Swedish electronic computer BESK at the Swedish Board for Computing 
Machinery (MNA) in Stockholm, where the program is available under the de- 
signation IM1/0732. The program is, of course, not restricted to the use of SCF 
functions, but can equally well be applied to Slater functions. In the latter 
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case, however, a more convenient program, extremely simple in use, has been 
devised by Flodmark [10] and recently presented in this journal. It is available 
at MNA under the notation S3. 


APPENDIX I 


Additional data regarding the program 


A single precision (9 place) floating point interpretive routine named Flinta 
was utilized in the program. This diminished the space available in the com- 
puter, causing the following limitations to be imposed on the number of the 
quantities used: 


Maximum number of k,, y, and 6, values =48 different of each 
c, values =18 different 
A, (y) and B,(6) integrals =10 of each (0<y, » <9) 


The limitations set to w and », essentially introduced to avoid the use of se- 
decimal numbers, give rise to the restrictions on n4+n, given in the text. 


Computation of A,,(y) and B, (6) 


A, (y) was computed with the aid of the following relationship 


A, (y) == Ag-1(y) + Ao (y). 


a3 RS 


In computing the B,(6) integral, two types of formulae were necessary as men- 
tioned in the text 


Ae 6. = 9 
B,(6) =" By 1 (8) +! a : |5|>4 
1 ° 5! 
a a . 16/<4 
Bel?) Si eoa are Tay aIpEne foes eke 
5 5 5° | 
‘ . & dd, |6|<4 
ar\g) a Re IEEICES fo 13 


As can be seen the value |5| has been used instead of the correct one |6/| 
as a criterion for the choice of formula to be employed in computing B, (6). 
This approximation, however, was found to be quite satisfactory. 


APPENDIX II 


Referring to the example in the text the contribution to the overlap inte- 


gral can be written 
1 


Sy= [ vapode=22 {| o(69) Xba exp (—8y-n0s) dE dy 
g=1 n=-1 4 


= Sy (A, (y) By (0)). 567 
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From Fig. 1, z=24— 40. Introducing z4= 14 cos 34 = $0 (1 + €n) giving z= 3 o&n, the 
corresponding moment integral is expressed by 


oo Jd 
My | yaryode— S22 | [ en-aen 
€=1 n=-1 
x > ky, exp (— Eyn— 1 dx) dE dy =5 * Sy (Aus (y)* By +1 (6). 
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